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Abstract 

Recently the Belle collaboration discovered two resonances, ^^(lOGlO) and ^^(lOGSO), that lie 
very close to the BB* and B*B* thresholds, respectively. It is natural to suppose that these are 
molecular states of bottom and anti-bottom mesons. Under this assumption, we introduce an 
effective field theory for the ^^(10610) and Z(,(10650), as well as similar unobserved states that 
are expected on the basis of heavy quark spin symmetry. The molecules are assumed to arise from 
short-range interactions that respect heavy quark spin symmetry. We use the theory to calculate 
line shapes in the vicinity of B^*'>B^*'> thresholds as well as two-body decay rates of the new bottom 
meson bound states. We derive new heavy quark spin symmetry predictions for the parameters 
appearing in the line shapes as well as the total and partial widths of the states. 
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Recently the BELLE collaboration observed two resonances, Z'b(10610) and Zb(10650), 
in the decays T(55') — t- T(r?,5')7r+7r^ {n = 1,2, or 3) and T{5S) — ?■ hh{'mP)'n-^7i^ (m = 1 or 
2) pp. The ^(,(10610) and ^^(10650) (which we will refer to as and below) have widths 
of about 15 MeV, and their masses lie a few MeV above the BB* and B*B* thresholds, 
respectively. However, an analysis in Ref. [2] concludes that an interpretation of the states 
as bound states lying below the B^*^B* threshold is still consistent with the available data 
on the decays T{5S) /ifjTT+TT . This conclusion depends on using line shapes 

for the that account for the coupling of the Z^'^ to the nearby B^*^B* thresholds rather 
than the Breit-Wigner form that was used in the experimental analysis. The experimental 
analysis favors the quantum numbers I^{J^) = l^(l^) for the Zf, and states. Arguments 
based on heavy quark spin symmetry [HI 0] indicate that there should be similar states called 
Wbo and W^q with quantum numbers I'^{J^) = l^(O^), as well as possibly Wbi and Wb2 
with quantum numbers l^(l^) and 1^(2+), respectively. 

In this paper, we will assume that these states are weakly bound molecules of heavy 
mesons. This is the approach adopted in Refs. [SH7|. For alternative interpretations of these 
states as tetraquarks, see Refs. P4T0]. Ref. [25] uses existence of the X(3872) and arguments 
based on heavy quark symmetry to argue that molecular states in the bottom sector must 
exist. Assuming the states are weakly bound molecules means they can be studied using a 
low energy effective field theory (EFT) that consists of nonrelativistic kinetic terms for the 
mesons and contact interactions whose coefficients are tuned to produce the bound states 
with energies close to threshold. A theory of this kind called XEFT has been developed for 
X{3872), which is thought to be a shallow 5-wave bound state of D^D*'^ + D^D*^ [nHll]. 
The purpose of this paper is to construct the analogous theory for the isovector Z^, Z^, and 
W^] states. This theory is similar in structure to the pionless effective theory used for very 
low energy nuclear physics [151 IIS]- K will be used to derive line shapes for the resonances 
that are valid near the relevant thresholds as well as calculate the two-body decays 

of the resonances. The predicted line shapes and decay rates incorporate the constraints 
imposed by heavy quark symmetry. New predictions for the parameters in the line shapes 
and the total and partial widths of these states are obtained. Experimental tests of these 
predictions should aid in interpreting the newly discovered Z^ and Z'^ states and searching 
for their partners. 

The EFT of this paper can be applied when the relative momentum of the B^*^B^*^ 
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mesons is much smaller than the pion mass. In two-nucleon scattering, a box diagram with 
two-pions is suppressed relative to one-pion exchange by Q/Ajy^ where Q p uIt, and 
= g\MN/{8nf^) ^ 1/(300 MeV), where Mn is the nucleon mass, gA is the nucleon 
axial coupling, and is the pion decay constant ^17]. Perturbative treatment of pions fails 
rather badly in two-nucleon systems [H] when p > 171.,^. For pion exchanges between a B* 
and B^*^ meson, the expansion parameter is Q/Abb, where 1/Abb = q'^Mb/ {^t^ fl)- Here 
Mb is the B meson mass and the axial coupling of heavy mesons, g, is between 0.5 and 0.7, 
which yields 160 MeV < A^b < 320 MeV. Since A^^ ^ Ajvtv we expect that perturbative 
treatment of pions will fail in the B meson sector for p ~ m^, but a pionless effective theory 
should work for p <^ rriT^. 

The Lagrangian we will use to describe low energy B^*^B^*^ scattering is 

^ = Tr[//i Udo + ^ J H,] + ^TrlHl a^] (1) 

V / ba 

V / ab 

- ^TT[Hlri,Hl,Hbr^,Hy] - ^TT[Hlri,a^Hi,HbnUH,] . 



Here a and b are SU(2) isospin indices, isospin matrices are normalized as t^t^^ = S^^ , 
traces are over spin indices which are not explicit, and Ha{Ha) is the heavy meson (heavy 
anti-meson) superfield. In terms of components. Ha = Pa + V^cr^ and Ha = Pa ~ KjO"', 
where Pa {Pa) and V^* {Va) are the pseudoscalar and vector B (B) mesons, respectively. The 
transformation properties of the fields under heavy quark spin and other symmetries are 
given in Ref. [12]. In Eq. ([T]), the mass M in the kinetic terms is the spin-averaged B 
meson mass, M = {3Mb* + Mb)/ 4: = 5314 MeV. The first three terms are the leading 
heavy hadron chiral perturbation theory Lagrangian of Refs. [T9H2T], written in the two- 
component notation of Ref. |22]- The next three terms are the Lagrangian for anti-heavy 
mesons. The kinetic terms must be promoted to leading order to prevent pinch singularities 
m 

scattering, as is conventional in nonrelativistic theory. The heavy mesons and 
anti-heavy mesons interact via the remaining terms in the Lagrangian which are contact 
interactions that mediate S-wave heavy meson scattering. Contact interactions of this type 
first were written down in Ref. [23], where the operators considered are proportional to 



Tr [if^ifj Tr[-ff^^b] and Tr [if^iJaCr*] TT[Hla^Hh]. It is easy to see that these operators can 
be written in terms of the single trace operators given above using Fierz transformations. 
In the case of scattering, one can classify states in terms of the total spin of the 

heavy quark and antiquark {Sqq), and the total angular momentum (Sqq) and isospin (/) 
of the light degrees of freedom. Since the light degrees of freedom in the B and B meson 
are isodoublets and have spin-1/2, the possible states of the light degrees of freedom are: i) 
/ = Sgg = 0, ii) / = and Sgg = 1, iii) 1=1 and Sgg = 0, or iv) / = Sgg = 1. The four 
operators in Eq. ([T]) mediate S-wave scattering in each of these channels and the notation 
for the coefficients is that the operator with coefficient Cjs mediates scattering in the isospin 
/ and spin S = Sgg channel. In the heavy quark limit scattering should be independent of 

Ref. [3] classified possible bound states of B and B mesons and concluded that there 
should be at least four and maybe six such isotriplet states. The wavefunctions of these 
states in terms of their Sqq and Sgg quantum numbers are derived in Ref. (See also 



Eq. (25) below.) Wbi and Wb2 are pure states with Sqq = Sgg = 1. The remaining states are 



mixtures of Sgg = or 1 and Sqq = or 1. If the mechanism that leads to shallow bound 
states operates in the Sgg = 1 channel or both Sgg = and Sgg = 1 channels, one expects 
to find all six shallow bound states. If the mechanism operates only in the Sgg = channel, 
then only four shallow bound states {Zj,, Z^, Wbo, and VF^'q) are expected. We will primarily 
focus on the former case, then comment on the latter at the end of the paper. 

The interpolating fields for these states are given by (we will drop the subscript b in what 
follows) : 

Z^' = ^{VlriPb-PaT^X) (2) 

U/AX _ Xyi Ayj 

where efj is a basis for symmetric traceless polarization vectors normalized as e^e<|^' = 5^^' , 
and a and b label fiavor antifundamental and fundamental indices, respectively. It is also 



possible to define isoscalar interpolating fields that are obtained from those in Eq. ^ by 
dropping the index A and replacing — )• 5ab/y/'^- We will focus on isovector states in what 
follows, the generalization to isoscalars is straightforward. It is enlightening to rewrite the 
contact interactions in terms of these interpolating fields. For the isovector fields these are 



^contact = -2 Cn (w.^'^W,^' + J2 (3) 



2 ^ M v^(^ii - Cio) Cio + 3Cn J \ 

\ C*ii — Cio C*io + Cii I \ Z ^ 
= -2 {wo+W^_, + ^Zf + W^'^W^' + ^ W^/.Vs^ ) (4) 

where in the last line the interactions are diagonalized by defining the fields Wqj^ = ^Wq"^ + 
^W^, = ^W^^-lW^ and = j^{Z^±Z^'). The Lagrangian for isoscalar terms 
is obtained by dropping the superscripts A and replacing Cu — )■ Cqj, z = or 1. Though we 
have diagonalized the interactions in Eq. Q, we will not work in this basis because the Z"^ 
and Z'^ are split by the hyperfine splitting, A = 46 MeV, and the and Wq^ are split 
by 2A = 92 MeV. 

It is straightforward to calculate the T- matrix for B^*^B^*^ scattering in these channels. 
For the channel we find 

where T,^*§*{E) is computed from a one loop diagram containing nonrelativistic B* and B* 
propagators of total energy E and is given by: 

M 



S 



B*B 



(E) = ^ (a - VM{2A -E)- te) . (6) 



The energy E is measured with respect to the BB threshold. The linear divergence in 
Eq. ^ is cancelled by the coupling constant 



so the T-matrix is given by 

~ 'M + ^//(2A -E)-Te ■ 

The T-matrix has a bound state pole at = 2A — jf-^^/M for 711 > 0. A similar calculation 
for the channel yields 

4vr 1 

Tw^ = TJ = , (9) 

M + ^M(A -E)-ie ^ ' 

which has a bound state pole at = A — If shallow bound states and W2 exist, 

heavy quark symmetry predicts their binding energies to be the same. On the other hand, 
if 7ii < 0, then there are no shallow bound states. In this case, heavy quark symmetry 
predicts the 5'-wave scattering length for B meson scattering in these channels to be the 
same. 

For and Z'^ states we must solve a coupled channel problem. The T-matrix is given 

by 

Tz' = -Ct - ^ziE) , (10) 
where Cz and Sz(-E) are matrices given by 

Cz = I + -^A, (11) 

Cii — Cio CiQ + Cii J 

^z{E) = I '''^'^^''^ ' ] (12) 

J^BB'iE) J 

M ( A - ^M{2A -E)-ie 
4^ \^ A - v^M(A -E)-ie 

^BB^E) = ^ (a - VAf(A -E)- le) . (13) 



The cutoff dependence in the T-matrix can be completely cancelled if the coupling Cio has 
the same form as Cu in Eq. ([T]), i.e., if Cio = C'io(A) = 27r/ (M(— A + 710)). For the T-matrix 
we find 



(14) 
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with the components given by 



2. _ ^ 47r -7+ + v/M(A - - i6 

M - ^M(A - - «e)(7+ - ^M{2A - E) - it) - 7^ ^ ' 

Tz'z = Tzz< = — — ^ (16) 

M (7+ - yM(A -E)- ie)(7+ - yM(2A - £;) - ie) - H ^ ' 

^ ^ 47r -7+ + VM(2A - i?)-i6 

M (7+ - ^M(A-E) -?e)(7+ - v/M(2A - - ^e) - 7^ ' ^ ^ 

where 7± = (711 ± 7io)/2. 

Solving the analogous problem in the Wq and Wq channels, we obtain 

™' M {^w _ ^_ME - ie){^Y - V^(2A - E) - ie) - (7^)2 ^ ^ 

47r 75 
W = W = M (^^ - v/=M^^)(7f ' - yM(2A - £;) - ie) - (7?^ ^^^^ 

y ^ 47r -77' + VmK^E)^e 

M {^w _ ^_ME - ie)(7f' - ^M{2i^ - E) - le) - (7^)2 ' ^ ' 

where 7^ - (710 + 37n)/4, 7^' = (8710 + 7ii)/4, and 7^ = ^(711 - 7io)/4 = ^87-/2. 

These amphtudes are valid near the relevant threshold. Specifically, if we take the molec- 
ular state to have a binding energy of 7^/M, then for 7 <^ m,r we can expand the amplitudes 
in powers of 7/-\/MA, where -\/MA — 494 MeV. For example, the state in the vicinity 
of the B*B* threshold will have energy E = 2A — 7^/M. Expanding the denominator in 
the expression for Tz'z' to leading order in 7, we find 

Tz,z> = — ^ , (21) 

M -^^ + 7 + 0(7VVAM) 

The pole in the amphtude is at 7 ~ 7+ or £^ ~ 2 A — 7+/M, corresponding to a Z' 
mass mz' = 2Mb* — iX/M- For the other states the binding energies are ^\/M, where 
A=Z, Z',Wo,Wq, Wi or W2-, and are given by 

Iz = Iz' = 7+ (22) 

IWi = Iw-i = 7ii 

710 + 3711 7z + 7vFi 

IWo = ^ = ^ 

3710 + 7ii 37z - 7wi 
= i = 2 ■ 

These relations between the binding momenta are consequences of heavy quark symmetry. 



We can incorporate the effects of decays of these resonances on the hne shapes by exphcitly 
violating unitarity. If the decays to other states also respect heavy quark spin symmetry, 
then we expect that incorporating these decays will just give imaginary components in the 
couplings of Eq. (jl]). So we promote Cqo, Cqi, Ciq and Cu to complex values, which also 
means 700, 701, 7io and 711 are complex. For each of them, we can write 775 = —l/ajs + 
iTis/2, where ais is the scattering length and Tjs is the total width of the bound state in 



the IS channel. The relations in Eq. (22) will still hold since they are a consequence of 
heavy quark spin symmetry and now give relationships among their imaginary components, 
i.e., the total widths. One prediction 

Til = nw^] = T[W2] = ^nWo] - Inw;,] , (23) 

was first derived in Ref. [Ij. In addition we also find that 

(24) 



r+ = -(En + Eio) = m = T[Z'] = -{T[Wo] + E^]) . 

The relation T[Z] = E [Z'] was first derived in Ref. |3j , while the last equality of Eq. ( p4| ) is 
new. One could derive this result, as well as similar predictions for partial decay rates, from 
the following decomposition of the wavefunction of the molecular states in terms of their 
components of definite Sqq (g) Sgg 



W2 
Wi 



Wo: 
Z' : 

Z : 



^, 

2 

2 
1 

1 
71 



^11 



.7=2 



J=l 



(25) 



QQ 



^QQ 



^qq + 2 ^QQ 



j=o 2 
L.- - i 1 



'-qq 



^QQ 



J=o 

Oqq 



•-qq 



Uq + 



^^QQ 

1 



Jqq 



V2 







qq ■ 



The molecular states inherit their widths from those of their constituent states with definite 
Sqq ® Sgg. Since the same constituent state appears in multiple molecules, by restricting to 
decays which are sensitive only to one choice of Sqq (g) Sgg, one can arrive at relations among 
molecular decays. Eor Sqq <^ Sgg = 1 ® 1, the result is Eq. (23). To derive the prediction in 



Eq. (24), consider the two S^ 

1 



2 

QQ 



\s 



QQ 



1,5^ 



0) 



spin configurations of the two heavy quarks: 

1 



V2 



(in) + i;t)), |5, 



QQ 



o,si 



0) 



V2 



(Iti)-lit)). (26) 



If Sq is the operator that measures the magnetic quantum number of the heavy quark only, 
then it is clear that 

2^^|1,0) = |0,0), 2^^|0,0) = |1,0). (27) 

Now let Mg^s' denote the interpolating fields with Sqq = s and Sgg = s', with all other 
indices labeling other quantum numbers suppressed for compactness. Then, it follows from 
the above that 

[2^^, Ml,] = Ml,, [2S'q, Ml,] = Ml,. (28) 

In what follows, h stands for any bottomonium state with allowed quantum numbers and 

Sqq = 0, and i is any allowed configuration of light hadrons. We define = 2SQ\h), which 
means /i is a bottomonium state with Sqq = 1. Then for the matrix element mediating the 
transition W, hi, we find 

{hi\W;,^Q) = ^{h£\Ml,\0) (29) 



{hi\[2S^Q,Ml,]\0) 



2 

{he\2S^QMl,\0) 



2 

Ch(''='k\Mlf='^\o) 



2 

3 



^ -y|(/i(5^=o)£|Z't(5^=o)|0) 

Rotational symmetry can be used to extend this result to other values of S^ and moreover 
implies that 

^)^(S-=m)^|^t(5^=m')|o) oc 5""""'. (30) 
Applying the result for h = r]b, which means /i = T, it follows that 

l-M(W^o ^ = lxlll\MZ^ ri)\'. (31) 

spins 

A similar analysis can be performed starting with a Wq instead. We find 

nWo^Vb^] ■■ nK^Vbi] ■ nZ^Ti] ■. rlZ'^re] = 1 : ^ : l : l. (32) 



This result is valid in the extreme heavy quark limit, in which T, rjb, Z^''^ and Wj"* are all 
degenerate. ^ In reality decay rates will also depend on the available phase space, which 
can be sensitive to the hyperfine splittings in each of the multiplets. For example, in the 
decays to single pions calculated below, the rates are multiplied by a pref actor of A;^ or 



and these kinematic prefactors introduce large corrections to Eq. (32) when the splittings 
between multiplets are not significantly larger than the hyperfine splittings of the multiplets. 
For decays to P-wave bottomonia, a similar derivation shows that 

^ Xbii] : rK ^ Xbii] : r[Z ^ hi] : T[Z' ^ hi] = ^ : ^ : 1 : 1 , (33) 



in the heavy quark limit. Eqs. (32) and (33) together imply the relationships among total 
decay rates in Eq. ( 24 ) assuming decays to quarkonia dominate the decays of the molecular 
states. 

One can explicitly check the claimed relationships among decay rates by computing the 
rates to final states with one pion using HHxPT. To obtain these rates, we will add to the 
Lagrangian of Eq. ([T]) the following interactions 

>ChhxPT = 9 Tr[Hl a'H,]Ai, - g Tv[HIH, a^] (34) 

+ ]^g^r,n^A^iHaH,]Al, + ^(?x,n tr[Tt + h.c. 

1 ■ - i 

The first line above gives the -kB^*^ and -kB^*'^ interactions. The next two lines are the 
interactions of bottomonium states with the B-mesons. Heavy quark spin symmetry groups 
the S- and P-wave bottomonium states into the multiplets 

T„ = aiT\nS) + 7]k{nS), (35) 
Xn = ^i{x'i{nP) + l^e^^^xTMP) + ^5\,,{nP)) + h{nP). 

^ The degeneracy of Z^'^ and W'^'^ is not a consequence of heavy quark spin symmetry, e.g., the binding 



energies predicted from Eq. (22) are not the same. However, these corrections to the masses of these 
resonances are 0(7^/Af ) and are expected to be a few MeV or less, and hence small compared to mass 
differences inherited from the hyperfine splittings of their constituent mesons. This lack of degeneracy is 
because the Z*^') and W'^'^ are linear combinations of members of different heavy quark spin multiplets, 
cf. Eq. (25). The multiplets of heavy quark spin symmetry are the (.^+, W^) and [Z^, W^) defined in 



the last line of Eq. ([3| and in terms of these states heavy quark spin symmetry predicts T\W — rii,l] 

T[z- re] and r[w+ xbi] = r[z+ ^ hbi]. 
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Note that the Lagrangian is heavy quark spin symmetric except for the hyperfine sphtting 



terms proportional to A, so the symmetry is restored in the hmit A = 0. From Eq. (34), 
one can compute the decay rates using the methods given in Refs. [TTHT4] . The decay rates 
to S-wave bottomonia and a single pion are listed below. We will drop the quantum number 
n, labeling the radial excitation level, but it is important to keep in mind that the coupling 
constants (^ttT, Q-wxi 9t and will be different for distinct multiplets. The decay rates to 
S-wave bottomonia and a single pion are 



r[Z ttT] 

r[z' ttT] 



8-KmwJ^ 
nirKEl 



E^{E^ + A) 
kl 1 



X Ci 

A m2 



(36) 



E^\ 3E^-A 



X Co 



2ggr§^ii 

kl 



AE^ 



3 ^2 - A2 
1 A 



2Am2 2 

+ t: 



A 



ggr 



The decays to P-wave bottomonia and a single pion are 



2 2 

+ - 



ggr 



kl A 
EIE^-A\ J 



AJ 

2i 



XO3 
X O4. 



T[Z ^ vr/ife] 
V[Z' -Kh] 
V[Wi ^ Txxm] 

T[Wi ^ 7rx62] 
T[W2 ^ vrxw] 
r[1^2 ^ vrx62] 



^h, kl 
12'Kmzfl 
^h, kl 

^XbO ^l 

^mx,,kl 
727rmvyJ2 

^xtikl 

mxt2kl 
8-Kmw2fn 



g^x + 



'^ggx 



E^ + A 

'^ggx 



E^-Ai 
1 



g-nx + 
g^x + 9gx' 
g^x + ggx 
g-Kx + '^ggx 
g^x + 



X C»i 

2 

X O2 

1 

E^ + A 
1 1 

3 1 



(37) 



A 



AE^- A 

E„ 



' X O4 

1 1 

\E^^-A)\ 



gnx + 
g-Kx + 



'^gg^ 



2 



E^ + A 

'^ggx 

E^-Ai 



X Ofi 



g^x + 



'^ggx 

Ett 



X a. 



Decays not listed here can only proceed through higher-derivative interactions which are 
suppressed compared to those listed. For compactness, we have defined the following non- 
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perturbative factors: 



O2 

06 



PariP,\W^)? 



(38) 



V2 



15 



Each of these has a prefactor of 1/3 to prevent overcounting when summing over the isospin 
states A, and an additional prefactor of 1/(2 J + 1) for a molecular state with spin J to 
prevent overcounting when summing over the spin polarization. In the limit of exact heavy 



quark spin symmetry, using arguments like that given in Eq. (29) one can show that the 
Cj's are all equal. 



The ratios of the decay rates in Eq. (36) and Eq. (37) match the predicted results from 



Eq. (32) and Eq. (33), respectively, when A = and E^^ is the same for all decays, which will 
be the case when the heavy quark spin symmetry multiplets (T^, r/fe), (x6j, hi,), and the Z^''^ 
and Wq^ are degenerate. Our explicit calculations of the decay rates allow us to incorporate 
important corrections to heavy quark spin symmetry predictions that come from phase space 
and kinematic factors. In HHxPT there are two mechanisms that contribute to the decay 
of the bound states. There is a short- distance process, mediated by the contact interactions 
g-j^r and g.,^^, in which the transition to the final state quarkonium and pion at a 

point. If this process dominates, the predicted ratios of rates are the heavy quark symmetry 



predictions in Eqs. (32) and Eqs. (33) weighted by factors of -E^/Ctt and fcjj, respectively. 
There is also a long-distance process 

^ B^*^B^*h followed by coalesence of the 
B and B meson into the final state quarkonium through the couplings gr and g^^. These 
processes lead to a more complicated dependence on the pion energy. When these processes 
dominate, the dependence is well approximated by kl/E^ since in all cases E.,,':^ /S.. 

For decays to S-wave bottomonium, the two processes appear at the same order in the 
expansion. If contact interactions dominate, which is obtained when the dimensionless ratio 



12 



15 








10 








5 




..J 


i 











-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 

FIG. 1: Ratio of T[Wo r]biT] (dashed), r[W^ ^ r/fovr] (dotted) and T[Z ^ Tvr] (solid) to 
r [Z' — )• Tvr] as functions of g gx /g-rrr 

r (MeV) 




FIG. 2: Partial decay rates T[Wo ??b(35)7r] (dotted), r[VF^ ^ ??fe(3S)^] (solid), T[Z T{3S)tt] 
(dashed) and T[Z' — T(35')7r] (dashed-dotted) as functions of g gr / g-KT for one possible choice of 
currently undetermined parameters {O g'^y = 10^^). 

of couplings, At = ggr/dnr ^ 1, then the ratios of partial decay rates are predicted to be 

r[Wo 7r7]b{3S)] : r[W;^ irVbi^S)] : r[Z 7rT(35)] : r[Z' 7rT(35)] 

= 0.26 : 2.0 : 0.62 : 1 (At = 0) , (39) 

where all partial decay rates have been normalized to the rate for r[Z' — ?■ 7rT(3S')]. In the 
opposite limit, 

r[Wo 7rr]b{3S)] : r[Wl^ 7^Vb{3S)] : r[Z 7rT(3S)] : r[Z' ^ 7rT(3S)] 

= 0.12 : 2.1 : 0.41 : 1 (|At| = oo) . (40) 

In computing these ratios, we have assumed Oi = O2 = ■ ■ ■ = Oq holds without significant 
corrections from symmetry violating terms. The values used for the masses of the bottomo- 
nium decay products, mr{3s) = 10355 MeV and ^^^^(35) = 10328 MeV, were determined in 
Ref. using a relativistic quark model. 
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The results do not differ greatly because ratios of the kinematic factors and k^/E^ are 
not that different. For intermediate values of Ax the results can depend rather dramatically 
on Ax when Ax ~ 0.6, as seen in Fig. [1} However, this wild variation occurs because of a 
cancellation between the two processes mentioned above. For these values of Ax, all four 
partial rates are highly suppressed, as shown in Fig. |2} In this figure, we have chosen 
Oi = (lOOMeV)^ and gl^ = l.OGeV"^ and calculated the rates as a function of Ax. Our 
choices of Oi and g'^j- are based on naive dimensional analysis and are only intended to be 
estimates that should be accurate within a factor of 10. For this particular choice of Oi and 
g^Y^ the branching fraction for Z — t- Ttt is less than 1%, for any value of Ax close enough to 



0.6 that the ratios deviate significantly from those given in Eqs. (39) and (40). Though this 
branching fraction has not been measured, the observation of this decay leads us to expect 
that the parameter Ax does not take on values where such cancellations suppress the decay 
rates. Therefore, we expect that experimental measurement of the ratios will yield results 



close to those in Eq. (39) or Eq. (40). 



For decays to P-wave bottomonium, the processes mediated by contact interactions are 
suppressed in the power counting. The relative importance of leading order to contact 
interaction mediated processes is controlled by the dimensionful parameter g^^/g^ which we 
expect to be ~ 1 GeV~^. In the limit g-Kx/dx = we find 

^ nXbi{2P)] : r[iy^ ^ 7CXbi{2P)] : T[Z ^ nh{2P)] : r[Z' ^ nh{2P)] 

= 0.72 : 0.57 : 0.66 : 1 {g^J g^ = OGeV-^) , (41) 

and 

r[lVi ^ 7rxw(2P)] : T[W2 ^ 7rxw(2P)] : ^T[Wo ^ rfXbi{2P)] - ^ r[iy^ ^ 7rxw(2P)] 
= 0.81 : 1 : 0.43 {g^Jg^ = GeV-^) . (42) 

The masses used to compute these ratios are m^^^{2P) = 10233 MeV, my^^^(2p) = 10255 MeV, 
"^Xb2(2P) = 10269 MeV and mh^^^p^ = 10261 MeV. The first three of these came from the 
Particle Data Group and the last is computed in Ref. [21]. For g^^/ g^ ^ 1 GeV~^ these 
ratios change by only a few percent. For g^^/ g^ — (300GeV)~^, there is a cancellation 
between contact and leading order diagrams which suppresses the total rates and which 
leads to modifications of the ratios similar to what was discussed above in the decays to 
5- wave bottomonium. Again the observation of Z' — )■ hb{2P)n disfavors such a suppression 
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of the decay rates. Such a cancellation between leading order and next-to-leading order 
contributions is also inconsistent with the power counting of the theory. 

Throughout this paper, it is assumed that all six isovector molecular states exist and are 
loosely bound. As stated earlier, it is also possible that binding occurs in the Sqq = channel 
only so that only the Wq ^ and Z^'^ bound states exist. In this case, the treatment of this 
paper is still applicable. One simply takes the binding momentum 711 = l/an to be small 
but negative. In this case, resumming both the Ciq and Cn interactions remains necessary. 
There is also the possibility that in these channels there are no shallow bound states or large 
scattering lengths. In this case, one can keep the summation in both channels as presented 
in this paper but tune 711 so there is no large scattering length. It should also be possible 
to sum only the strong, binding interaction and give a perturbative treatment to the weaker 
interaction. The T matrices computed using the latter approach should correspond to a 
power series expansion of those in this paper. We will save any investigation along these 
lines for future work. 

In this paper we introduced the Lagrangian describing heavy quark spin symmetric 5*- 
wave contact interactions among observed and hypothesized isovector B meson molecules. 
We derive the line shapes of these states in the vicinity of their respective thresholds, 
including coupled channel effects where mixture of states is possible. By doing so, we 
have arrived at relationships among the binding energies and decay rates of the molecular 
states. Some relationships among the widths of the Wj^] and Z^'^ states derived in this 



paper appear in Ref. [1], while the prediction in Eq. (24) is new. A confirmation of these 
predictions by explicit calculation of partial widths for strong two-body decays to S- and 
P-wave bottomonia using HHxPT was performed. This allowed us to compute corrections 
to the earlier predictions which arise from differences in the kinematics between the various 
processes. Tests of these predictions will aid in interpreting the new states. Future work 
could include an extension of these results to the isoscalar sector of B meson molecules, a 
detailed look at radiative decays, as well as a determination of currently unknown parameters 
Ax and gx/9nx using the angular distribution of decay products. 
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